arXiv:1504.03908v2 [hep-ph] 25 Jan 2016 


Dark Higgs Channel for FERMI GeV 7-ray Excess 

P. Ko and Yong Tang 

School of Physics, 

Korea Institute for Advanced Study, 

Seoul 130-722, South Korea 

Abstract 

Dark Higgs is very generic in dark matter models where DM is stabilized by some spontaneously 
broken dark gauge symmetries. Motivated by tentative GeV scale 7 -ray excess from the galactic 
center (GC), we investigate a scenario where a pair of dark matter X annihilates into a pair of 
dark Higgs H 2 , which subsequently decays into standard model particles through its mixing with 
SM Higgs boson. Besides the two-body decay of H 2 , we also include multibody decay channels 
of the dark Higgs. We find that the best-fit point is around Mx — 95.0GeV, Mh 2 — 86.7GeV, 
(cju) ~ 4.0 X 10“^®cm^/s and gives a p-value ~ 0.40. Implication of this result is described in the 
context of dark matter models with dark gauge symmetries. Since such a dark Higgs boson is very 
difficult to produce at colliders, indirect DM detections of cosmic 7 -rays could be an important 
probe of dark sectors, complementary to collider searches. 
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I. INTRODUCTION 


Firm evidences for dark matter (DM) come exclusively from the gravitational inter¬ 
action at the moment. A popular scenario in particle physics models, DM as weakly- 
interacting massive particles (WIMP), generally predicts that DM should have a mass be¬ 
tween (P(GeV) — (P(TeV), with a weak-scale annihilation cross section around 3x 10“^®cm^/s. 
If those annihilation hnal states go to standard model particles eventually, there might be 
notable excesses in cosmic rays and gamma ray searches. 

By analyzing Fermi-LAT's public data, several groups nmn have been claiming that 
there might be some excess in the gamma-ray signals from Galactic center inner Galaxy 
and even some Dwarf Galaxy [13]. The excess is at ~ O(GeV) energy scale and its 
morphology against the distance to galaxy center is consistent with signals from WIMP 
DM annihilation, which has motivated intense discussions about DM model-constructions 
and constraints [T3H72] . Besides DM interpretation, astrophysical origins of the excess have 
also been actively investigated |HQH90]. For instance. Refs. [HOl EH] showed that gamma-ray 
emission from unresolved millisecond pulsars is compatible with the excess and can account 
for at least part of the excess. And Ref. [89l |90] argued that the excess might be comprised 
entirely of point sources. Recently, Fermi-LAT collaboration published a paper |9T] on search 
for DM from Milky Way dwarf spheroidal galaxies and gave stringent constraints. However, 
since dwarf galaxies might have different DM density prohles from Milky Way halo, such 
constraints would be relaxed then. 

In this paper, we shall exclusively consider DM interpretations. As a hrst step, it is 
natural to investigate the GeV excess through annihilation channels that a pair of DM goes to 
two SM particles directly, such as qq, cc, bb, ti, gg, hh, WW, ZZ and also their different 
combinations with some branching fractions. After all, no new particle has been found yet 
at the LHG, except the Higgs boson. This method has provided valuable information for the 
favored DM mass and annihilation cross section ranges. Discussions has been extended to 
cascade two-body decay through new mediators, such as Z' and dark Higgs H 2 , which are 
ubiquitous in new physics beyond SM. In particular, light mediator {Mz\h 2 < IGeV) (SH] 
and heavy Z [2H] cases have been investigated thoroughly. 

This work is intended to investigate GeV scale gamma-ray excess in models where DM 
annihilates into a pair of heavy dark Higgs (> IGeV) which subsequently could decay into 
multi-body hnal states such as WW* or W*W*, where W* is a virtual W boson. The aim 
is to provide the ranges of the favored dark Higgs mass, DM mass and the annihilation 
cross section, which might be useful for particle physics model building, such as hidden 
sector DM models with local dark gauge symmetries. This work differs from many previous 
investigation in one essential aspect: we take into account consistently all possible decay 
modes for heavy dark Higgs, not restricted to its two-body decays. 

This paper is organized as follows. In Sec. we discuss the theoretical motivation and 
establish our formalism and notations. In Sec. in we present our numerical results on the 
best-ht parameters. Then in Sec. IV, we briehy discuss the implications for DM relic density 
and constraints from our results. Finally, we give a summary. 


^ Recently, Fermi-LAT also released an paper with some excess [Hj. 
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II. FORMALISM 


We shall consider the following annihilation channel for self-conjngate DM X, 


X + X ^ H 2 + H 2 , followed by 

H 2 SM + SM{+SM). 


Here H 2 denotes the dark Higgs, distinguishing it from the SM-like Higgs Hi with — 125 
GeV. In Ref. [37] , the present authors showed that dark Higgs is very generic in dark matter 
models with dark gauge symmetries, and the GC 7 -ray excess can be easily accommodated 
with the above mechanism (see also Refs. [711[33] for related works). 

H 2 can decay into SM particles through its small mixing with Hi. The mixture between 
H 2 and H 2 can be easily achieved in particle physics model building. 

For example, a real scalar DM X and a complex scalar $ (dark Higgs that breaks local 
dark gauge symmetries) can have the following interactions. 


C 3 


( 2 . 1 ) 


where H is the SM Higgs doublet. After gauge (or even possible global) symmetry breaking. 


0 , 


Vh + h 

TT 


-)■ 


v^ + (p 


( 2 . 2 ) 


where Vh and are the vacuum expectation values, two neutral scalars h and 0 will mix 
with each other through the Higgs portal coupling \^h, resulting in two mass eigenstates 
Hi and H 2 with 

\H 2 J \ sincK cos a J\vJ 
in terms of the mixing angle a. 

The above Lagrangian is just one example of many DM models with Higgs portal. One 
can also consider the case with a real scalar 0 , 


C D -AiXV - _ x^(t)H^H - \i(t)‘^H^H. (2.4) 


Again after the electroweak symmetry breaking. 


0 , 


Vh + h 

h and 0 are mixed. We can also consider a model with fermionic X with 

C D -A 1 X 75 X 0 - X2(j)H^H - Xscp'^H^H. 


(2.5) 


( 2 . 6 ) 


All the above models can easily evade current experimental bounds when a is very small 
(see discussion in Ref. [92] for example). If one assumes all other possible new particles are 
heavy, DM X will dominantly annihilate into H 2 ^s. 

To be as general as possible, we shall just work with the effective operator, X^i7| 
(X 75 XWI for fermionic X), and consider the annihilation process in Fig. assuming that 
other particles in the dark sector are all heavy enough. 
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FIG. 1: Feynman diagram due to the effective operator for fermionic X or 

X^X^H 2 for vector X). The actual annihilation process may occur through s or t channel, and 
contact interaction. Details in the gray bubble depend on various ultraviolet completions. The 
produced H 2 S can have two-, three- or even four-body decay channels. 


The produced 7 ^ 2 ’s could be either relativistic or non-relativistic for Mh^ Mx or 
Mh 2 — Mx, respectively. Different kinematics might lead to signihcant differences in the 
gamma-ray spectra. Moreover, depending on the mass of H 2 , Mh 2 , H 2 can dominantly decay 
into 2 or 3 standard model particles. In our numerical calculation, we use PYTHIA-6.4 [HS] 
to simulate and tabulate dN^ / dE^ for the interesting ranges of Mx and Mh 2 ■ In particular, 
we focus on Mx > 5 GeV and Mh 2 > 1 GeV. 

The general differential flux of the gamma-ray from the annihilation of self-conjugate DM 
is given by 


dE^dVL 


1 ^^7 

8 vr^ Mly^ dE^ 


drp^ {r' (r, 0 )), 


(2.7) 


where {o'v)l^^ is the velocity-averaged annihilation cross section for the annihilation channel 
/, dN^/dE^ is prompt gamma-ray spectrum, r' = \/rQ + — 2rQr cos 9, r is the distance to 

earth from the DM annihilation point, Tq ~ 8.5kpc for solar system and 6 is the observation 
angle between the line-of-sight and the center of Milky Way. An extra factor 1/2 needs to 
be included for non-self-conjugate DM, such as complex scalars or Dirac fermions. In our 
considered case, we have only one annihilation channel, X + X ^ H 2 + H 2 . 

For DM density distribution, we use the following generalized NFW profile [M] . 


with parameters Vc — 
not stated otherwise. 


P = P© 
20 kpc and p© ~ 


r®i 

7 

'1 + rQ/rc 

- r . 


1 + r/vc 


0.4GeV/cm^. We shall adopt the index 7 


( 2 . 8 ) 

1.26 if 


III. NUMERICAL ANALYSIS 


We first show three cases for the gamma-ray spectrum in Fig. The vertical axis marks 
the conventional 


MN . 1 /■ , 

_ _ rpZ _ / _ /o i\ 

dE ~ XVL J^^dE^dVt' 

where AD indicates the region of interest. The 24 data points we used to compare with are 
from Ref. m, denoted as GGW hereafter. 
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FIG. 2: Three illustrative cases for gamma-ray spectra in contrast with CCW data points m- All 
masses are in GeV unit and av with cm^/s. Line shape around E ~ Mji/ 2/2 is due to decay modes, 
H 2 77 , Z'y. 


As we can see, different parameter sets can give different spectrnm shape, especially in 
the high energy regime. When the branching ratios of H 2 —>■ 77 , Zj are increasing, we can 
see the gamma lines more easily aronnd E ~ Mh 2 /‘^- Since the annihilation cross section is 
at order of 10 “^®cm^/s and the branching ratios of H 2 —)■ 77 , Z'y are around 0 . 2 % at most, 
the considered parameters are still consistent with constraint from gamma-line searches. 
We now use the function and hnd its minimum to hnd out the best £t: 

x" (Ma-, (av)) = /•) (ft - fi ). (3.2) 

where /Xj and fi are the predicted and measured fluxes in the i-th energy bin respectively, and 
E is the 24 x 24 covariance matrix. We take the numerical values for fi and S from CCW m- 
Minimizing the against fi with respect to Mx, Mh 2 and (av) gives the best-fit points, 
and then two-dimensional Icr, 2a and 3a contours are defined at Ay^ = y^ — Xmm = 2.3, 
6.2 and 11 . 8 , respectively. 

Fig. [^is our main result. In the left panel, Mx, and (av) are freely varied, so that 
the total degree of freedom (d.o.f.) is 21. The red dot represents the best-fit point with 

Mx — OS.OGeV, — 86.7GeV, (av) ~ 4.0 x 10“^®cm^/s, (3.3) 

gives Xmin — 22.0, with the corresponding p-value equal to 0.40. 

We also notice that there are two separate regimes, one in the low mass region and the 
other in high mass region. The higher mass region is basically aligned with Mh 2 — Mx since 
otherwise a highly-boosted H 2 would give a harder gamma-ray spectrum. In this region. 
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FIG. 3: The regions inside solid(black), dashed(blue) and long-dashed(red) contours correspond 
to Icr, 2a and 3cr, respectively. The red dots inside la contours are the best-fit points. In the 
left panel, we vary freely Mx, Mh 2 and {av). While in the right panel, we fix the mass of H 2 , 
Mh 2 — Mx- 


H 2 mostly decays into bb. As one increases the mass of H 2 , H 2 —)■ H 2 —t 

H 2 —t 77 and H 2 —)■ 'jZ become more and more important, and all of them give harder 
gamma-ray spectra either due to the leptonic hnal states or the gamma lines. Eventually, 
increases significantly when Mh 2 > ISOGeV. 

In the low mass region, the contours are scattered but centered around Mh 2 — lOGeV 
and such a light H 2 most likely decays into bb, cc and . Since cc and t^t~ would give 
harder spectra than bb does, we would need a lower Mx to fit the data, which is exactly 
what we see in Fig. (dotted curve). Increasing the branching ratios of cc and would 
require a even lower Mx- 

In the right panel, we show a special case in which Mh 2 — Mx, so that the produced H 2 S 
are non-relativistic. In such a case, the d.o.f. is now 22. The best-fit parameters are 


Mx — Mh 2 — 97.1GeV, {av) ~ 4.2 x 10 ^^cm'^/s. 


(3.4) 


which gives ximn — 22.5 and the p-value equal to 0.43. An interesting thing is that Mx — 
Mh 2 — 125GeV also give a good-fit. This point is equivalent to the channel that DM X 
annihilates into SM Higgs, which has been already found in previous study [sniisa. 

In the left panel of Fig. we £x the mass of dark Higgs to the best-point value, Mh 2 = 
86.7, and vary Mx and {av). We show la, 2a and 3a contours in terms of solid(black), 
dashed(blue) and long-dashed (red) curves, respectively. To compare with bb channel, we 
also present 3a region in the right panel of Fig. The best-fit point is around 


Mx — 49.4GeV, {av) ~ 1.75 x 10 ^®cm^/s, 
which gives Xmin — 24.4 and a p-value, 0.34. 


(3.5) 
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FIG. 4: Regions inside solid(black), dashed(blue) and long-dashed (red) contours correspond to Icr, 
2cr and 2>a, respectively. The red dots inside la contours are the best-fit points. In the left panel, 
we hx Mh 2 = 86.7GeV, but vary freely Mx and {av). For comparison, in the right panel we 
consider XX ^ bb channel. 


Channels 

Best-fit parameters 

Xmin/d-O.f. 

p-value 

AA ^ H 2 H 2 
(with Mh 2 / Mx) 

Mx — 95.0GeV, Mh 2 — 86.7GeV 
{av) ~ 4.0 X 10“^®cm^/s 

22.0/21 

0.40 

XX H 2 H 2 
(with Mh 2 = Mx) 

Mx — 97.1GeV 
{av) ~ 4.2 X 10“^®cm^/s 

22.5/22 

0.43 

XX HiHi 
(with Mhi = 125GeV) 

Mx — 125GeV 
{av) ~ 5.5 X 10“^®cm^/s 

24.8/22 

0.30 

AA ^ bb 

Mx — 49.4GeV 
{av) ~ 1.75 X 10“^®cm^/s 

24.4/22 

0.34 


TABLE I: Summary table for the best fits with three different assumptions. 


For completeness and comparison, we also show similar plots in Fig. based on the 
analysis of Inner Galaxy from Ref. ra. All parameters are the same as previous except 
7 = 1.18 for the DM profile as used Ref. [lO]. As shown in Fig. the regions inside contours 
are much smaller than those in Figs. and This is due to the fact that uncertainties 
in Ref. [10] are purely statistical. All the best-fit points give ^ 44, and the minimal 
is reached when Mx — 87.0GeV, Mh 2 — 82.1GeV and {av) ~ 4.7 x 10“^®cm^/s with 
Xmin/d-o.f — 42.6/(22 — 3) which corresponds to a p-value, 1.5 x 10“^. It seems that this 
is not a good fit. However, according to Ref. [10], any value of < 50 should be taken 
as a reasonable “good fit”, given the large systematic uncertainties associated with the 
background templates choice. Nevertheless, we can easily see that the favored regions are 
consistent with those in Figs. and [^ within 2-3cr. 
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FIG. 5: Similar plots as Figs. and based on results from Ref. m- Regions inside solid(black), 
dashed(blue) and long-dashed(red) contours correspond to la, 2a and 3a, respectively. The red 
dots inside Icr contours are the best-fit points. 


IV. RELIC ABUNDANCE AND PHENOMENOLOGY 


So far, we have not discussed any actual concrete DM models, how the correct relic abun¬ 
dance can be achieved and how one can test or constrain this scenario further. Generally, 
such topics are highly model dependent and the conclusions would differ significantly from 
one case to another. Still, we can give some general implications from the results obtained 
earlier in this letter. 


From Eqs. (3.3) and (3.4), we notice that the best-fit annihilation cross section for dark 
Higgs channel is bigger than the canonical value for thermal DM A, 2 ~ 3 x 10 ^®cm^/s. 









although still consistent within 2-3a range. Larger (av) would mean a smaller relic abun¬ 
dance for X. However, the correct relic density can still be reached if we extend the above 
minimal DM model setup. For example, we may introduce another heavy DM component 
Y and the total energy fraction of Dy -|- Dx is just the observed Ddm- Suppose Y froze out 
in the early universe but had decayed into X pairs with a lifetime shorter than the age of 
Universe. This can be easily achieved if we introduce the following interactions: 

6C = yiY^X^ + y2MxYX^. 


When y 2 is small enough, say less than 10“^^, Y would decay into X pair after BBN epoch or 
the freeze-out of X’s. Before its decay, U’s relic abundance was determined by yi through the 
efficient annihilation process, Y + Y X + X. Although j /2 is small, it is still technically 
natural or natural by ’t Hooft’s naturalness argument because if j /2 is zero, we have an 
additional Z 2 symmetry for particle Y. _ 

For bb channel, the best-fit cross section, Eq. (3.5), is a little smaller than the canon¬ 
ical value. This is not a problem at this stage since the right panel of Fig. shows that 
canonical value of thermal cross section is within l-2a range. Even if it is a problem, we 
can still imagine that in some concrete models, there could exist some other annihilation 
channels whose cross sections have velocity dependence or p-wave suppression, (av) oc 
For instance, if fermionic DM X has interaction, XXcfZ, then X -|- X —)■ 20 is much more 
suppressed nowadays than at the freeze-out time. Then the sum of cross sections to 00 and 
bb has a canonical value at the freeze-out time but only bb channel is important now for its 
indirect detections. 

Here we briefly discuss the related possible particle phenomenology in an explicit UV 
complete model. We shall note phenomenology is highly model-dependent and we refer to 
our previous works [SQ EZl EHl ES for some detailed examples. Our discussion here is just 
focused on the following Lagrangian with f/x(l) gauge symmetry: 

1 


£ = £sm - ($t<F - u|/2)' 


- vl/2) {H^H - vl/2) 


(4.1) 


where in the SM Lagrangian the Higgs potential term is \h [H^^H — and the covariant 

derivative is defined as — igxX^)^. 

A nonzero vacuum expectation value vq> breaks U{l)x spontaneously. Afterwards, 
gets mass Mx = gxVq> and the dark Higgs field 0 will mix with the SM Higgs field h through 
the Higgs portal term, resulting in two mass eigenstates. Hi and H 2 . The mixing angle a 


in the matrix Eq. (2.3) is determined by 


sin2Q; = 


2\}j^vhv^ 


m 


Hi 


m 


(4.2) 


H, 


In this model, gauge boson is the dark matter because its only couplings, X^X^iLf 2 
X^X^iLi_ 2 ; have accidental Z 2 symmetry which ensure the stability. 

The mixing angle is constrained by Higgs signal strength and invisible decay at the 

LHC unmnsi, 


sin^ a < 0.2, and \h<s> ^ 10 ^ for < 62.5GeV. 
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FIG. 6: Feynman diagrams. The last one can be neglected for small A$. 


Due to the allowed small mixing, the dominant annihilation processes for small scalar 
self-coupling are shown in Fig. their thermal cross section (an) to determine the relic 
density is given by 




gtcos'tt 8(M^,-4Mi) 16 + 6MP 

144,rM2 [ 


(4.3) 


For the above formula, we can determine the required gauge coupling gx for correct DM 
density. For instance, if Mh 2 ^ Mx, (uv) ~ 3 x 10“^®cm^/s gives 


0.2 r Mx \ 

cos a ylOOGeV/ 


While if Mh 2 — Mx, {<Jv) ~ 3 x 10 ^®cm^/s would give 


0.24 / Mx \ 
cos a ylOOGeV ) 


The scalar mixing also lead to possible signal for direct detection of DM. The spin- 
independent scattering cross section for exchanging scalar mediators is calculated as 


SI _ / I 1 y 

” MlJ 




2.2 X lO-^^cm^ 


/ 9X^a^a 

V 10-2 


2 


/75GeV 

V Mh2 



(4.4) 


which can easily satisfy LUX |106j limit for small mixing angle sin a < 0.05. This also means 
that direct detection can place a much stronger limit on sin a for light dark higgs. Since the 
production cross section of this second higgs goes like sin^ a, it would be very challenge for 
probe it at the LHG. 


V. SUMMARY 

In the letter, we have explored a possibility that the GeV scale y-ray excess from the 
galactic center is due to DM pair annihilation into a pair of dark Higgs, followed by the dark 
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Higgs decay into the SM particles through its small mixing with the SM Higgs boson. Includ¬ 
ing the correlations among different parameters and varying Mx and Mh 2 independently, 
we ffnd that the best fit is obtained if 

Mx — 95.0GeV, Mh 2 — 86.7GeV, {av) ~ 4.0 x 10“^®cm^/s 

If we impose Mx — Mh 2 , we get a similar result (see Table I). This information could be 
important inputs in dark matter models with dark Higgs boson. 

At this stage we cannot make any strong statement about the existence of dark Higgs 
with mass close to the DM mass ~ 95GeV. However dark Higgs is very generic in DM 
models where DM is stabilized by some spontaneously broken local (or even global) dark 
gauge symmetries [silEaETiiniiEsHina]. Since the dark Higgs boson is a SM singlet scalar, 
it is very difficult to find it at colliders although there are some interesting constraints from 
the LHG data on the SM Higgs signal strengths |in4] . It is simply more difficult to produce 
them at colliders when the mixing angle is small. It is very amusing to notice that indirect 
DM detection experiments can be more sensitive to such a dark Higgs than the collider 
experiments. Gompared with more popular dark photon scenario, it is more natural to have 
flavor dependent couplings of dark Higgs boson to the SM fermions, since its couplings are 
basically the same as those of the SM Higgs boson modulo the mixing angle effect. This 
fact makes much easier DM model building [211122] . It remains to be seen whether this fit 
survives in the future data sets, and if there would be any indication of such a dark Higgs 
from the future collider experiments. 
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